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Abstract 

We express the unitary time evolution in non-relativistic regularized 
quantum electrodynamics at zero and positive temperature by a Feynman 
integral denned in terms of a complex Brownian motion. An average over 
the quantum electromagnetic field determines the form of the quantum 
mechanics in an environment of a quantum black body radiation. In 
this non-perturbative formulation of quantum electrodynamics we prove 
the existence of the classical limit Ti — > O.We estimate an error to some 
approximations commonly applied in quantum radiation theory. 
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1 Introduction 



Quantum electrodynamics (non-relativistic as well as relativistic) has been for- 
mulated in terms of a path integral by Feynman himself |l|. Starting from 
the Feynman integral the conventional perturbation series in powers of the fine 
structure constant a — e 2 (hc)^ 1 is derived. It is a common attitude among the 
physicists to treat the Feynman integral just as a tool for a generation of the 
perturbation series. However, it is well known that the imaginary time (Eu- 
clidean) version of the Feynman integral has a mathematical meaning beyond 
the perturbation theory. In fact, the Euclidean two-dimensional quantum elec- 
trodynamics (QED) has been constructed by means of the path integral(Q ||). 
The path integral for Euclidean regularized QED at finite temperature in any 
dimension has been discussed in ref.[||. Although, in principle, the Euclidean 
correlation functions can be continued analytically to the Minkowski space || , 
in practice, this is an inconvenient approach to dynamics. 

Perturbative QED agrees perfectly well with experiments. Its non-perturbative 
effects have not been discovered yet. However, the photons form an inevitable 
environment for all events in nature at least in the form of the heat or (much 
weaker) cosmic microwave background. The permanent influence of radiation 
may have a profound effect on mesoscopic phenomena. A lot of theoretical 
studies has been done investigating the effect of an infinite environment on 
a finite microscopic system (see e.g. || Q] Q). However, the mathematical 
model of the heat bath could not be identified with any physical environment. 
The quantum electromagnetic field constitutes such a physical environment. 
There are suggestions that Leggett's model approximates the effect of the 
quantum electromagnetic field. However, we think that this is quite a rough ap- 
proximation completely out of control. The approximations commonly applied 



[g_0j |12j (I5| |lj [[L3| ||f§ [|17| neglect the spatial dependence of the electromagnetic 
field. The neglect of the spatial dependence is hard to justify in QED where the 
electromagnetic field has singular vacuum correlation functions. 

In this paper we apply to QED our earlier mathematical formulation of the 
Feynman integral (T^-Jl^l (see also |2(J |^l|). Such a formulation makes it pos- 
sible to discuss the dynamics directly in the real time. Working with physical 
time is indispensable in problems concerning finite time dynamics e.g. decoher- 
ence |22| and radiation reaction. We consider the Dirac model of a finite 
number of particles interacting by means of the Coulomb forces and additionally 
with the quantum electromagnetic field. This is the first approximation to the 
relativistic quantum field theory of interacting quantum electromagnetic field 
(in the radiation gauge) and matter fields. As pointed out already by Feynman 
[0 a minor modification of the path integral defines the complete relativistic 
model. For this reason the non-relativistic theory has a similar perturbative 
expansion in powers of the fine structure constant a. At the moment we are 
mainly interested in an application of the quantum electromagnetic field as a 
model of an environment to the conventional non-relativistic quantum mechan- 
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ics. For this reason we include a discussion of mixed states corresponding to a 
finite temperature. The main virtue of a rigorous formulation of the Feynman 
integral is an opportunity to estimate an error to formal calculations. In this 
paper we discuss the radiation damping and quasiclassical expansion. In the per- 
turbative approach to QED we encounter a serious difficulty with the classical 
limit (see |p3| ), because a tends to infinity when Ti — > 0. So, a non-perturbative 
formulation of this problem is required. 

The paper is organized as follows. In sec. 2 we define the model in terms of the 
Hamiltonian in a finite mode approximation for the quantum electromagnetic 
field. In sec. 3 we derive the Feynman formula for the time evolution operator in 
this model. The finite temperature QED is expressed by the functional integral 
in sec. 4. In sec. 5 we integrate over the quantum electromagnetic field in order 
to describe the quantum mechanical particles in equilibrium with a radiation 
background. We consider the limit of an infinite volume and infinite number of 
modes in sec. 6. We define the Feynman kernel for the evolution operator in sec. 7. 
We prove a semiclassical asymptotics for the Feynman kernel in sec. 8. In sec. 9 
we discuss the approximation of a slowly varying quantum electromagnetic field 
(the dipole approximation) |lo| JI^] |Q . Finally, sec. 10 contains a summary 
and perspectives for further research. 

We put more effort on the precise description of quantum dynamics than on 
the details of the proofs. This is so, because the method of the probabilistic 
representation of the quantum dynamics has been described in our earlier papers 
[|l8]-(l^] together with detailed proofs. Our approach in this paper starting with 
a finite mode approximation is a repetition of methods of our earlier work. 

2 The Hamiltonian 

We consider the electromagnetic field in a finite volume V.We expand the vector 
potential A in eigenmodes g(k,x) 

w(fc) 2 g(fc,x) = c 2 Ag(fc,x) 

In this formula uj(k) denotes the eigenfrequencies and k numerates the eigen- 
modes (so k is is a vector if V is a parallelepiped ; this will be tacitly assumed 
further on, but we write some formulae in a way suggesting cavities of a general 
shape). We assume that g(k,x) are real (see e.g. refs.|f25f |2(| for such a choice 
of eigenmodes) normalized as follows 

/ g(k,x)g(k',x)dx = 5{k,k')V 
Jv 
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The (5-function is equal to 1 for the same g and equal to zero if the eigenmodcs 
are different. The expansion in eigenmodes takes the form 



In eq.(l) £ is the polarization transversal to k 

!/=l,2 

a(fc, i/) and a(k, v) + are the amplitudes fulfilling the (Poisson brackets) com- 
mutation relations for creation and annihilation operators. f(k) is a formfactor 
regularizing the electromagnetic field. Such a formfactor has been introduced 
already by Pauli and Fierz p7| (we choose f(k) real; with no regularization 

/(*0 = 1). 

The free Hamiltonian is 



H R = V ^(fc)a(fc, i/)+a(fc, i/) = V -4ttP (fc, + ^/j,(k)w(k) 2 y(k, vf (2) 
where we defined 



with 

^)-^ 

In this way we have expressed the radiation Hamiltonian as a sum of the Hamil- 
tonians for harmonic oscillators. The Hamiltonian describing an interaction of 
the r-th particle with an external potential u r , the interaction with photons (the 
coupling 2t , where e r denotes an electric charge and c the light velocity) and the 
Coulomb interaction among themselves (these are the longitundinal degrees of 
freedom of the electromagnetic field, hence the form of the regularized Coulomb 
interaction is a consequence of eq.(l))reads 

Hp = (p r + — A(x r )) 2 +^e r e s 0(x r ,x s ) + ^u r (x r ) (3) 

r r,s r 

where 

0(x r ,x s ) = 4ir^2f(k) 2 c 2 uj(ky 2 g(k,ic r )g(k,x s ) 

k 

When the number of modes is finite then the Hamiltonian Hp + Hp is defined as 
a self-adjoint operator in a tensor product of a finite number of L 2 (dy) Hubert 
spaces. 
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3 Probabilistic representation of quantum dy- 
namics 

If we restrict ourselves to a finite number of modes and quantize the Hamiltonian 
H = Hr + Hp canonically then we obtain a second order (in x r and y(k, v)) dif- 
ferential operator. A probabilistic representation of the unitary group generated 
by H has been discussed in our earlier papers fig)] . 

Let us consider first a single harmonic oscillator. Denote its Hamiltonian by 
Hk- The ground state \k of Hk is 

x k (y) = eM-^§^) (4) 

Let us define the Brownian motion as the Gaussian process with the covariance 

E[b t b s ] — min(t, s) 
Let q(t, y) be the solution of the stochastic equation 

dq = —iuj(k)qdt + Xa^db (5) 
with the initial condition y. Here 

*->+<> 

and 

/ h 

Ok 



fl(k) 

Then, the Schrodinger time evolution of the state ip = Xk® can be expressed 
in the form 

My) = (exp(-iH k ^)(y) = X k(y)EMq(t, y))] (6) 

The free electromagnetic field is just a collection of harmonic oscillators. Let x 
be the ground state of Hr i.e.Hpx — O.Let again rp — x'f-Then 

A(A) = (e X p(-iH R ^)(A) = x (A)EMA(t, A))] (7) 



where 



and 



A = ^ y/4:TT^(k)cf(k)£(k, v)y{k, v)g{k, x) 



A(t, A) = £ V^Mk)cf(k)£{k, u)q{t; y, k, v)g{k, x) (8) 
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where 



q{t;y ) k,v)—exp{—iuj(k)t)y(k,u) + \ak / exp(—iuj(k)(t — T))db(r;k,v) 

Jo 

is the solution of eq.(5) with the initial condition y(k, v) (now each Brownian 
motion has an index numerating the eigenmodes of the electromagnetic field; 
the Brownian motions corresponding to different indices are independent). 

In order to derive a probabilistic representation of the Schrodinger evolution 
for Hp + Hp we need first some additional cutoffs. So, let us assume at the 
beginning that the number of modes is finite. Next, we need to bound the be- 
haviour of A(x) as well as </>(x r , x s ) for complex x.In ref. JlJ] we have bounded 
A and <fi m the Feynman integral multiplying these functions by a Gaussian 
damping factor depending on the Brownian motion. We introduce here a sim- 
pler regularization particularly useful in QED. We regularize the modes g(k,x) 
as follows 

g(k, x) -» g s (k, x) = exp(-5 2 \x\ 8 )g(k, x) (9) 

A with the gs modes replacing g does not fulfill the sourceless Maxwell equa- 
tions. For this reason the replacement g —> gg is treated as a regularization which 
subsequently is to be removed. The regularization (9) can be related to the one 
of our earlier paper Jllfl by means of the Laplace transform 

exp(-6 2 \x\ 8 ) = / dsh(s)exp(-sVd\x\ 2 ) 



where h(s) is determined by the inverse Laplace transform (we made use of 
(iy) 4 = y 4 ) 

1 f°° 

h(s) = — dyexp(-y 4 + iys) (10) 

27r J-oo 

It follows from eq.(10) that h is a real, square integrable and bounded function 
on R. 

We obtain now the following Feynman formula 
Theorem 1 

Consider the Hamiltonian Hp + Hp (eqs.(2)-(3)), where the number of modes 
is finite and g(k, x) are replaced by gs(k, x)( this Hamiltonian is denoted Hs). Assume 
that gs(k,K) can be analytically extended to the region Vc 

V c = (z e C 3 : z = x + (1 + i)y, xeVjei? 3 ) 

and remain bounded analytic functions in Vc ■ Assume that the external poten- 
tials u r are also analytic and bounded in Vc-Let us consider the initial state of 
the form x$ where x is the ground state of Hp and $ is square integrable. Then, 
the unitary semigroup Iff (t > 0) generating the Schrodinger time evolution has 



G 



for a sufficiently small time t the probabilistic representation 
(f/ t <V)(A,x) = (eip(-iff a *)^)(A,x) = X (A) 



exp[iYLr TTc Jo As(t, x r + Xa r b r (T))Xa r db r (T) 
ex p{ ~ H Sr, s e »- e s Jo <M X >- + Acr r b r (r),x s + Xa s b s (T))dr ' ' ' ' 

-fE r Jo M*r + Xa r b r (r))dT) <t>(A 5 {t, A), x + Xah{t)) 



Remarks 

f .The index 5 at the potentials A and <p means that in their definitions we 
made the replacement g — > g§. 

2. It does not matter whether the stochastic integral in eq.(10) is in the 
Stratonovitch or Ito sense Q because both integrals are equal for transverse 
vector potentials. 

3. If V is a parallelepiped then cos(kx) and sm(kx) constitute the eigen- 
modes, hence gs(k,x) are analytic and bounded in Vc- 

4. The wave function in the theorem describes the quantum electromagnetic 
field interacting with the n-particle system. We denoted the coordinates in the 
argument of ip and $ collectively by x and b. 

5. We may write Ut — (U±) n , where the natural number n is chosen large 
enough . In such a case we can obtain a probabilistic representation for arbi- 
trarily large time through a repetition of the formula (11). 

6. We consider here only t > the evolution for t < can be obtained by 
time reflection as discussed in ref. jl8| . 

Proof: the proof is based on general methods of the construction of semi- 
groups [^9|.The factor in the expectation value in eq.(ll) is an example of a 
non-anticipating multiplicative functional. Such a multiplicative functional de- 
fines a semigroup. There remains to prove that the multiplicative functional is 
integrable (so that the expression (11) is finite) and that the generator of the 
semigroup is — j-Hg.The integrability results from our assumptions implying 
that 4>s and As are bounded. Then, for each As 



exp^z J As{t, x + Xab(T))Xadb(T) 



is integrable for any complex z on the basis of an estimate proved in ref. |]30| and 
discussed in our earlier paper Q. There remains the Gaussian integral over 
As- Owing to the (5-regularization the covariance of As is a bounded continuous 
function. Hence, the ,4,5-integral is bounded by an exponential of a continuous 
bilinear form in the Brownian motion b (this bilinear form is discussed in more 
detail in sec. 5). The expectation value of such an exponential is finite for a 
sufficiently small time (see |}2| Jl8) ) . We can refer now to the Fubini theorem in 
order to perform first the ^4,5-integral and to conclude that the expectation value 
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in eq.(ll) is finite. Finally, the computation of the generator involves the limit 
t — ► 0. In such a case the calculations are reduced to an elementary stochastic 
calculus pHI which does not depend on the fact whether the potentials are real 
or complex. Let us explain in more detail only the quantum electromagnetic 



part which constitutes a new ingredient in comparison to ref. 1 18 . In the differ- 
entiation of Ag(t,x) over t we expand As in q(t] k, i/).Then, the differentiation 
of q is determined by eq.(5).As a result of the q(t; k, v) differentiation we obtain 
the photon part of H 

x ~ 1Hrx = - g imm^M + M%(fc ' v) oah)) 

The remaining terms result from differentiation in time in the same way as in 
the conventional form of the Feynman-Kac formula p3] . 

Let us mention the imaginary time (Euclidean) version of eqs.(5) and (11) 
(the Feynman-Kac formula). The Feynman-Kac formula can be applied to more 
singular potentials. However, there are stronger conditions on the possible de- 
crease of the scalar potential to minus infinity. The imaginary time version can 
be useful in a study of the spectrum, but is unsuitable for an investigation of 
dynamics. So, the imaginary time version of eq.(5) reads 

dq(t; k, v) — ~u>(k)q(t; k, v)dt + Ukdb(t; k, v) (12) 

whereas the imaginary time version of eq.(ll) takes the form 

(exp(-Hj-)i;) (A,x) = 



X (A)E 



expi 



exp(^ - i J2 r ,s e r e sJo 0(x r + cr r b,.(r),x s +a s b s (r))dT 
~k Jo u r(x r + <J r h r (T))dT) <S>{A(t, A), x + ab{t)) 



(13) 



where A(t, x) is expressed by q(t; k, v) in the same way as in the real time case 
(eq.(8)).Eq.(13) holds true under the assumption that the number of modes is 
finite and g(k, x) are bounded functions on R 3 (in such case no <5-regularization 
is needed). It is easy to take the limit of an infinite number of modes if the 
formfactors f(k) decay sufficiently fast for large k (e.g.|/(fc)| < const\k\~ 2 ). 

4 Probabilistic description of quantum dynam- 
ics at finite temperature 

A mixed state can be described as a random pure state. Let us begin with one 
degree of freedom. Assume that the (classical) probability p n of a state |n > is 
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given. We introduce a random variable w such that 

E[w n w m ] = 5 mn n\p n 
Then, the density matrix can be expressed in the form 

P = E[\ X XX\] (14) 

where 

oo 1 

| x >=£ w >> (15) 



Let us consider the quantum harmonic oscillator with the frequency u)(k) at 
temperature T . If \n > is the n-th excited energy eigenstate of the oscillator 
then for the Gibbs distribution 

p n = (1 — exp(~[ihuj{k))) exp(—(3huj(k)n) 

where (3 — (KT)^ 1 and K is the Boltzman constant. It is easy to see that the 
expectation value (14) can be realized by the Gaussian distribution 

Tr(pF) = Tr (E[\ X >< x\F}) = (1 - exp(-0fiw(k))) Tr (^j dp k (w)F\ X >< x\ 

where dpk is the Gaussian probability measure 

dpk(w) = Tt~ 1 dwdwexp{f3hLo(k))exp{—exp(l3huj(k))\w\ 2 ) 

We can sum up the series (15) . In the coordinate representation we obtain an 
expression for the coherent state 

< y\Xk(z) >= { — )^exp{-\-q k y 2 - ^z 2 + r] k zy) 

where 

h 

Vk 



p(k)Lo(k) 
and we defined 

2 

z = \ j — w 



Vk 

\x(z) > is just a coherent state with the mean value of the position equal to 
Rez. Under the time evolution 

Xk(t, y) =< y\exp(-iH k ^)\x(z) >=< y\x(exp(-iu(k)t)z) > (16) 
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According to our general scheme |28j the time evolution of the wave function 
|x > (eq.(15)) determines the correlation functions in the mixed state p. In 
more detail, we define first the stochastic process (0 < r < t) 

dq(r) = — X (t - T, q)- 1 Vx(t - T, q)dr + Xadb(r) (17) 
m 

Then, the correlation functions of Heisenberg operators can be expressed by the 
solution q(t, y) of the stochastic equation with the initial condition y e.g. 



Tr(y(t)y(0)p) = E[J dy\ X t(y)\ 2 yq t (y)] (18) 

In the case of the coherent state Xk eq.(17) reads 

dq(r) = —iu(k)qdT + iuj(k)zexp(~iuj(k)(t — r))dr + \akdb(r) (19) 

Its solution entering the formula (18) is 

q(t; y, z) — exp(—ito(k)t)y + iexp(— iu(k)t)zsin(u>(k)t)+ 

+\a k Jj exp(-ioj(k)(t - T))db{r) ' 20 ' 

The correlation functions in eq.(18) can easily be computed from eqs.(16) and 
(20) 

Tr(y(t)y(0)pj = (l - exp(-0fiw(k))\ J dp k (w)dy\ < y\xk{exp{-iuj{k)t)w) > 
(exp(-iu(k)t)y 2 + sin(u(k)t)exp(-iu(k)i)wyj 

= 2^(k)u(k) (cth(±(3huj(k))cos{uj(k)t) - isin(u(k)t)j 

(21) 

According to the formulae (2) and (8) the free quantum electromagnetic field is 
a sum of oscillators. The density matrix is a product of the density matrices with 
different masses and frequencies. Hence, we obtain directly from eqs.(8) and (21) 
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7Y^A,-(i,x)A,(0,x>J =G /S (x,x';% = 
^Ek^^gik^giky) ^ (22) 
c£/i(i/37ic|k|)cos(c|k|i) — £sin(c|k|t) 



The expectation value in the ground state x of the free electromagnetic field in 
a cavity is a special case of eq.(22) corresponding to the limit (3 — » oo (t > 0) 



D jZ (x,x';t) =< X |^(t,x)^(0,x')|x >= 

f(fc) ; 

|k| 



^ J2k ^%9{k, x) ff (k, ^)exp{-ic\k\t) (23) 
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The expression for correlation functions of an n-particle system in equilibrium 
with the black body radiation at temperature T described by the Hamiltonian 
Hr + Hp (2)-(3) follows directly from our description (14) of a mixed state as 
a random pure state. In this way we obtain directly from Theorem 1 

Tr (u+^A, x)f/ t $ 2 (A, x)p) = / \[ k l/ dy(k, v)dx$i(A, x) 
E $2(A z s (t, A), x + \ab(t))exp(i £ r t So A I( T > x - + Aa r b r (T))Aa r db r (r) 
exp(^ - j- J2 r , s e r e s J * <f>s{Xr + Aoyb r (r),x s + \a s b s (T))d,T- 

-Ji Y.r So u r{*r + X(T r b r (T))dT 



(24) 

In eq.(24) A z means that in the expression (8) for A we insert the solution (20) 
depending on the coherent state random variable z^ as well as on the initial 
condition y(k, v). Then, in eq.(24) the average is taken over the Brownian mo- 
tion, the random variables Zk as well as over the initial values y(k, v) defining 
A. Note, that the left hand side of eq.(24) is simple if <&i = 1 whereas the right 
hand side is quite complicated. This is so, because the unitarity condition ap- 
plied at the left hand side when expressed by a path integral takes a complicated 
form (a condition of a time invariance of the expectation value is an analog of 
unitarity) . 



5 Quantum mechanics in QED environment 

In quantum mechanics we usually consider idealized systems consisting of a fi- 
nite number of particles. However, real systems are open and interact in an 
uncontrollable way with an infinite environment. It is believed that if the in- 
teraction with the environment is weak then it is negligible. A model suggested 
by Caldeira and Leggett JjJ indicates that this may be not the case concern- 
ing processes on the mesoscopic scale. In particular, the disappearance of the 
quantum interference in macroscopic systems may be a result of an interaction 
with an environment (besides the Caldeira and Leggett papers see also [[35| 
[|36|). However, the model of an environment discussed in the above mentioned 
papers does not come from any concrete physical environment ( the claims in 
refs.[p[]ll[| that such an environment can be considered as an approximation to 
QED is rather unconvincing). There appeared some suggestions |37) |H that 
QED (e.g. in the form of the black body radiation) constitutes the environment 
enforcing the appearance of classical properties of macroscopic systems. Un- 
fortunately, the quantum electrodynamics discussed in these papers has been 
considered in such a rough approximation that a conclusion concerning its ef- 
fect was rather unfounded. Another reason for a discussion of an effect of a 
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black body radiation in quantum mechanics is its possible substantial effect on 
Rydberg atoms. 

A restriction to observables associated with a system of particles means an 
average over all states of the quantum radiation field. Let us note first that the 
random electromagnetic field (8) consists of two parts:the ^-independent part 
(which we call A vac ) depending on the initial values y(k, v) 

A vac (t, x) = V^Kk)cf(k)£(k, v)exp(-iuj(k)t)y(k, v)g{k 1 x) (25) 
and the rest 

A Q = A - A vac 

An easy computation gives 
E[Af(r,x)A? = 

^ E fe ^f l9 (k, x) 5 (fc, x') (esp(-ic|k| |r' - r|) - exp(-ic\k\ (r' + r))j 
= D j7 (x,x'; |r' - r|) - £>Vf c (x,x';r' +r)) 

(26) 

where we divided the correlation function (26) into a time-translation invariant 
part and the rest. If the 5-regularization g — > gs is inserted in eq.(26) then the 
corresponding correlation function is denoted D s . 

We can calculate now transition amplitudes in quantum mechanics in a QED 
environment. Let x(A) be the Fock vacuum for the radiation field. Assume that 
$i and $2 depend only on the particle positions. Then the transition amplitudes 
can be expressed in the form 



$ 2 (x + Acrb(i)) 



expi - i Ers So So db r (r)db s (T')D s ( Xr + A<r r b r (r),x s + Aa s b s (r'); \t - r'\) 

expi - j- Ysr^ e r e s f Q 4>s(y~ r + A<7 r b r (r),x s + Xa s b s (T))dr- 

-|Er So u r{*r + Xa r b r (T))dT 



(27) 

where dhdhD means a sum over vector indices of D and vector indices of the 
Brownian motions of individual particles. 
Remarks: 

1. The conditions for applicability of the Fubini theorem are fulfilled in 
the regularized theory for a sufficiently small time (cp. Theorem 1). Double 
stochastic integrals are discussed in refs.|39 | 40[ . 

2. It is remarkable that the "vacuum part" in the exponential of the expec- 
tation value (27) is absent. The "vacuum part" of eq.(26) canceled with the 
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contribution from the expectation value of A vac and what remains is the time 
translation invariant part of the expression (26). 

The corresponding formula at finite temperature reads 



Tr^C/+$ 1 (x)L/ t $ 2 (x)pJ = J dx$ 1 (x) J B^ 2 (x + Aab(i)) 

exp( -iJ2 rs ^Sf^ f* f* db r (r)db s (r')G^(x r + Aa r b r (r),x s + A<r s b s (r'); \t - r'|)) 
expi - j- J2 r , s e r e s f Q ^(x r + Acr r b r (r), x s + ACT s b s (r))dr- 
_ t Jo "r(x r + Af7 r b r (r))dr 



where Go is the finite temperature correlation function (22) with g — > 35. 



(28) 
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We can see from eq.(27) that the quantum electromagnetic environment con- 
tributes an additional non-local interaction (described by D ) between particles 
of the system. One could say that the particles interact through an exchange 
of photons from the environment. As D (at zero temperature) is proportional 
to h the additional interaction is of the order of the fine structure constant 
a. From eq.(22) we can see that at large temperature the interaction becomes 
proportional to the temperature. However, because of the light velocity in 
the denominator, this interaction seems negligible at the room temperatures. 
Nevertheless, it should be realized that such a comparison of the strength of 
interactions makes sense only for interactions of the same type. The interaction 
induced by the QED environment is of different type than the Coulomb inter- 
action <fi and the external potential it. It depends on velocities and is non-local. 
Such interactions may be relevant on the large time and space scale. 



6 Removal of cutoffs and regular izat ions 

We restricted ourselves till now to a finite number of modes (an ultraviolet cut- 
off) in a finite volume V and we introduced additionally the 5-regularization 
which is demanded by the functional integral representation in the real time 
(it would be superfluous in the Feynman-Kac formula). Together with a finite 
number of modes we considered also the formfactor f(k) which can describe 
a finite distribution of the charge. We can allow the number of modes to in- 
crease to infinity if the formfactor decreases sufficiently fast (in fc).In such a 
case the formfactor fulfills the role of the ultraviolet cutoff. We restrict our- 
selves from now on to the case when V is a cube. Then, g(k,~x) consist of pairs 
(cos(kx), sin(kx)) (with a certain restriction on k ensuring that the basis is 
orthogonal, sec |||2j). Hence, 

g{k, x)g(fc, x') — > cos(kx)cos(kx') + sm(kx)sm(kx') = cos(kx — kx') (29) 
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It can be seen from eqs.(26)-(27) that in a cube (with periodic boundary condi- 
tions) the interaction is invariant under space and time translations. Moreover, 
if the temperature is equal to zero, the cube tends to R 3 and the ultraviolet 
regularization is removed (f(k) = 1) then the interaction becomes relativistic 
invariant ( in such a case D coincides with the conventional D + ).lt is well-known 
that in the relativistic theory D + is singular on the light cone. The origin of this 
sigularity can be seen from eqs.(26) and (29); the oscillating factor becomes a 
constant on the light cone, hence the sum over wave vectors tends to infinity. In 
such a case the double stochastic integral in eq.(27) is infinite. The ultraviolet 
problem is unavoidable in QED and presumably can only be solved perturba- 
tively through an expansion of the exponential in power series (which is a power 
series in a). We consider in this section the remaining regularizations. 
Theorem 2 

Let V be a cube . Assume that the number of modes is finite . Assume also 
that the external potentials u r are bounded on Vc- Then, the limit S — > of 
exp (—iH$lr) ip(A, x) exists for each A and x and sufficiently small t. The limit 
(5^0 defines a unitary evolution solving the Schrodinger equation 

with the initial condition ip. As a consequence the limits 5 — > of the transition 
amplitude (27) and the finite temperature correlation functions (28) also exist 
for a sufficiently small time. 

Proof: We apply the representation (9)-(10) to the 6- regularization. We 
expand the exponentials in eq.(ll) in power series in c.Thcn, it follows from the 
results of rcf.|Q (Theorem 3) and from the Lebesgue dominated convergence 
theorem (because As and 4>s are bounded) that for each A® the limit 8 — > is 
exchangeable with the sum. There remains the Gaussian integral over A® in the 
power series. From estimates of the combinatorial factor in front of the term 
(A®) n done in Jig] (based on ref.|3Q|) it follows that it is bounded by -^=. Hence, 

after the Gaussian integration of the term (A®) n this term becomes bounded 
by c(t) n where c(t) is proportional to t because the upper limit in the stochastic 
integration is equal to t. Hence, if t is sufficiently small then the power series is 
bounded by a convergent geometric series. The proof of the convergence of the 
power series for the transition amplitude (27) and the correlation functions (28) 
is similar. In the proof we replace the Gaussian integral with the covariance 
(26) by the Gaussian integral with the covariance (23). 
Corollary 3 

Assume that there exist a > and R(a) > such that 

|/(*)| < R(a)exp(-a\k\ 2 ) (30) 

i)Let N be the number of modes. Then, the limit N — > oo of the transition 
amplitude (27) and the finite temperature correlation functions (28) exists uni- 
formly in S for a sufficiently small time t. The limits N — > oo and S — > are 
exchangeable and can be taken simultaneously. 
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ii)If time t is sufficiently small then after N — > oo the limit V — > R 3 of the 
transition amplitude and the same limit of the finite temperature correlation 
functions can be performed according to the rule 



The limit V —> R 3 is uniform in 5. Hence, the limits V — > R 3 and S — > are 
exchangeable. 

Remark: We could define directly the evolution operator Ut (11) for an 
infinite number of modes and a sufficiently small time by the same argument 
which we apply in the proof below. In principle, an operator can be recovered 
from the transition amplitude. 

Proof:It is sufficient to note that with the bound on the formfactor f(k) the 
limits considered in the Corollary can be taken at each order of the perturbation 
series in e. This is so, because according to the results of ref.|l^] (proven under a 
condition equivalent to the bound on the formfactor assumed here) and Theorem 
2 the perturbation series is convergent. Then , the existence of the limit N — > oo 
is easy to check at each order of the perturbation expansion. The infinite volume 
limit and the convergence of the sum over k to an integral is a consequence of 
Riemann's definition of an integral. 

We approached QED starting from a well-defined Hamiltonian consisting of 
a finite number of modes of the quantum electromagnetic field in a finite volume. 
This is a model often studied in quantum optics. However, we could consider the 
final formula (27) where the limits 5 — > ,N — > oo and V — » R 3 are performed 
directly under the expectation value on the right hand side in the exponential 
factor. We would obtain an expression defined by the stochastic integrals of 
D and tj>. This could be considered as a definition of QED. In fact, if we first 
integrated over A in the (imaginary time) Feynman-Kac formula and afterwards 
performed in a formal way the analytic continuation on the path space, then we 
would arrive at the formula (27). In Theorem 2 and Corollary 3 we have shown 
that the limits on the right hand side exist for a sufficiently small time and 
define a limit of an evolution operator. If the limits are performed under the 
expectation value then the expectation value may be infinite or even if it is finite 
we may lose a connection with the Hamiltonian formalism. In the rest of this 
section we investigate the double stochastic integral of the formula (27) treating 
the right hand side as a definition of the Fcynman integral representation of the 
left hand side. We emphasize the following result of ref. ff(j| 

Lemma 4 

The sufficient conditions required for a definition of the double stochastic 
integral in eq. (27) as an almost surely finite random variable read 




E 



dr [ dr'\D (b(r), b(r'); r — r') | 2 < oo 



(31) 



Jo 
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and 

E ' 



dr\D(b(r)Mr)^ 



2 







< oo 



Applying Lemma 4 we can show 
Lemma 5 

Under the assumption (30) the double stochastic integral is defined as an 
almost surely finite random variable. 

Proof: Under our assumptions the correlation function D for equal times 
is a finite constant. There remains to prove that the integral (31) is finite. The 
expectation value (31)is equal to 

/„ dr/ r dT , (Px(PyS'kSk'{2Ti)- :i {T'{T - r'))"' 

-p(-^-^)(ikiik'ir i 

exp (-ic(\k\ + |k'|)(r - r') + Xa(ky - kx) + A<7(k'y - k'x)) 
We change variables (x, y) — > (x, z) where 



x y = \fr — t'-l 

After such a change of variables and a calculation of the x-integral the integrand 
becomes regular in r and r' whereas the integrability at large k and fc'is ensured 
by the assumption (30). 

Finally, we consider an explicit example 

/(fc) = eap(-||k|) (32) 

(where a > 0) which does not fulfill the assumption (30). Nevertheless, we can 
show 

Proposition 6 

Let us restrict ourselves in Theorem 1 to a single particle. 

i) If f(k) is defined by eq.(32) then the function D (the case of an infinite 
volume) has an analytic continuation to the region Vc(0) = (z G C 3 : z = 
(1 + i)y,y G R 3 ). Moreover, D is bounded in this region. 

ii) Assume that the external potential u is bounded by a quadratic form 
in the region Vc of Theorem 1. Then, the functional integral (27) with the 
formfactor (32) is finite for a sufficiently small time. 

Remark: The multiparticle case with the formfactor (32) is more involved, 
because D and 4> are meromorphic multivalued functions in Vc- 
Proof: A direct calculation from eq.(23) (see also eq.(26)) gives 

%(x,x'; T - T >) = M ( (a + lc \ T - T '\f + ( x _ x ') 2 )-i 
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For a single particle the function D in the double stochastic integral of eq.(27) 
takes the form 

D j7 (Aab(r),A<7b(r');T-r') = 5f t D = 

±5% (a 2 - c 2 (r - r'f + 2iac\r - r'\ + io\h{r) - b(r')) 2 )" 1 

8 tr is a bounded operator. So, it is sufficient to discuss the scalar function 
D behind S tr ( a more detailed demonstration would involve an application of 
the Ito formula ;we express the stochastic integral with respect to the longin- 
tudinal modes by an ordinary integral which would be bounded owing to the 
regularization). It is easy to see that D is bounded because, if 



then \D\ < (37ra 2 ) 1 because the real part of the denominator of D is different 
from zero. If 

\t-t'\c>- 

then \D\ < (3tto 2 )^ 1 , because the imaginary part of the denominator of D is 
different from zero. 

When D is a bounded continuous function and u is quadratically bounded 
then the integrand (27) can be bounded by the Gaussian integral. In such a 
case the well-known theorems (see e.g.[^2); b can be expressed by independent 
Gaussian random variables) show the integrability proving ii). 

Even if D is bounded it may be not easy to prove integrability in eq.(27) for 
arbitrarily large time. We have to show that the negative part of D decreases 
sufficiently fast when |t — T l \ grows. If this is not the case we have to continue the 
formula (27) beyond the time-interval where it loses its meaning. We encounter 
this problem already for a harmonic oscillator. A solution of the problem is 
well-known in this case . For an anharmonic oscillator it is discussed in a 
semiclassical approximation (Maslov indices). 



7 The Feynman kernel of the evolution operator 

We wish to derive an explicit Feynman integral formula for the evolution kernel 
of the operator (11). We mainly discuss in this section a finite number of modes. 
However, first we would like to explain that in principle there is no difficulty 
in defining an evolution kernel in infinite number of dimensions . There is no 
Lebesgue measure in infinite number of dimensions. For this reason we define 
a kernel K as a transformation of a set of measures Q into a set of functions 
L 2 (dfi) integrable with respect to a certain measure fi. In QED we choose for 
/i a measure defined on vector fields which gives a realization of the Fock space 
of the transversal electromagnetic free field i.e. /x is the Gaussian measure with 
the covariance 

MjK-A)" 1 
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If <&i and $2 belong to L 2 (d/i) then we define the kernel Kt of the unitary 
evolution Ut by the formula 

($ lf l7 t $ 2 ) = {$ U K t $ 2 dfl) EE / d 3 ^(C)$T(C,x) 

C, x; 0, B, x')$ 2 (B, ^^(B^V 

We have obtained a Feynman formula for the kernel of a quantum mechanical 
Hamiltonian evolution in the case without the quantum electromagnetic field in 
ref . jl9| . In order to extend this formula to a quantum electromagnetic field let 
us start with a definition of the Gaussian process < r < t 

Q*(r; k, v) = o*( T ; k, v) + Ao*7*(t; k, v) (33) 

where a* (we shall sometimes conceal some indices which are inessential at the 
moment) is the solution of the oscillator equation 

(|_2 +^(fe)VM=0 (34) 
with the boundary conditions a*(0) = X and a*(i) = Y", explicitly 

t f \ sin{uj{k){t - t)) sin(uj(k)T) 

a (T) = - — 7W~\ H 7 — f i \ \ 

sin(u){k)t) sin(uj{k)t) 

7* is the Gaussian process fulfilling the zero boundary conditions 7*(0) = 7*(i) = 
with mean zero and the covariance 

E[ 1 t {r) 1 t (T')]=M{r 1 T') (35) 

where 

C t M(t,t') = (-^2- oj(kf)M{ T y) = S(r - t') (36) 

i.e.CAi = 1.A4 is the Green function of the operator C with the Dirichlet 
boundary conditions on the interval [0,t].For t < iroj(k) the kernel M. is pos- 
itive definite, because for such a t the operator C is positive definite. This is 
so, because the lowest eigenvalue of — d 2 is (f) 2 - In such a case the random 
variable 7* is real. The explicit formula for M. reads 

M(t, r') - ^(r' - r) (sin(uj(k)(t - T '))co S {u{k){t - r))- 
—ctg(u)(k)t)sin((j(k)(t — T))sin(uj(k)(t — r') 

+ Zm ( T - T ') ( s ™M fe )(* ~ r))aw(w(k)(t - r'))- 
—ctg(u(k)t)sin(uj(k)(t — T))sin(u>(k)(t — r') 
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where 6 is the Heaviside step function. 

We can express the process 7* by the Brownian motion 

7 *(r) = sin(u>(k)(t- t)) J (sin{u(k)(t - u)jj 1 db(u) (37) 

We insert Q 1 (33) into the expansion (8) of the electromagnetic field A into 
eigenmodes. We obtain 

A(t, x; t, B, C) = A c1 {t, x; t, B, C) + AW*(r, x) (38) 

where A cl is the solution of the wave equation 

(^"A)^ = 

with the two-point boundary condition A cl (0, x; t, B, C) = B(x) and A cl (t, x; f , B, C) = 
C(x). 7v is a quantum fluctuation. Note that the two-point boundary value 
problem (34) has the unique solution only if t < for every k.In terms of the 
initial condition of eq.(34) B has the form 

B = y/*Mk)cf{k)£(k, v)X(k, u)g(k, x) 

We have a similar formula for C in terms of Y(k, v). The decomposition of 
the quantum field A (38) into a solution of a classical boundary value problem 
and a random fluctuation resembles Boyer 's stochastic electrodynamics pl| . 
However, our random fluctuation TL is multiplied by a complex number A. 

The construction of the random electromagnetic field (38) fulfilling the two- 
point boundary condition is a generalization of the construction of the Brownian 
bridge v |3^| . This is the Gaussian process defined on the interval [0,1] with 
zero two-point boundary condition and the covariance 

E[v j (j)v,(/)] - Sji(e(T - r')r'(l - r) + B(t' - r)r(l - r')) 

The bridge connecting the points x and x' which we apply in the Feynman 
formula has the form 

p(r;t,x,x') = x'^ + t —j-^x + \aV~tv(^) (39) 

We can now express the evolution kernel by the bridge 
Theorem 7 

i) Let < t < m i n k{^jh\) then the integral kernel of the unitary semigroup 
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Ut of Theorem 1 (under the assumptions of this theorem ) has the representation 



ii) assume that i) holds true then under the assumptions of Theorem 2 there 
exists the limit S — > of K s (t, C, x; 0,B, x') for each x,x',B and C if t is 
sufficiently small. 

iii) the limit V — > R 3 exists for a sufficiently small time if f(k) = for 
|k| > K with a certain ultraviolet cutoff k 

Remarks 

1. The allowed time in ii) and iii) may be smaller than in i) because we 
additionally require a convergence of the perturbation series in powers of the 
electric charge. 

2. The formula for arbitrary time can in principle be derived by a composition 
of kernels for a small time. 

Proof:In i) we can repeat the argument of ref. in order to represent the 
kernel of the evolution operator of Theorem 1 (eq.(ll)) in terms of the corre- 
sponding bridge joining the initial and the final points. Here, for the electro- 
magnetic field we apply the bridge (33) corresponding to the harmonic oscillator 
instead the w-bridge corresponding to a free particle used in our earlier paper 
| fl9l (the imaginary time version of the oscillator bridge is discussed in |Q).In 
part ii) we repeat the proof of Theorem 2 showing that the perturbation series 
is convergent uniformly in 6 and that the ^-regularization can be removed term 
by term. In part iii) we apply the same argument as in ii) of Corollary 3 to 
show that a sum over wave vectors k converges to an integral when V — ► O.This 
is possible only if the denominator in the bridge Q* (33) is different from zero. 
For this reason we assumed the sharp ultraviolet cutoff in iii). 

It seems that eq.(39) remains true for any t ^ ^xy- A rough argument is 
as follows: it could be proved in an independent way under our assumptions 
that the evolution kernel of Hr + Hp has a convergent expansion around the 
evolution kernel of Hr (see ref.[p4| for a closely related theorem). On the other 




(40) 
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hand eq.(40) defines such an expansion for any t ^ zj{k)' Hence, it is sufficient 
to compare the expansion coefficients. These, expansion coefficients must be 
the same because both kernels fulfil the same equation. Another argument 
which could be applied comes from the fact that an imaginary time version 
of eq.(40)(cp.eq.(13)) holds true for arbitrarily large time. For an operator 
bounded from below an analytic continuation in time exists (see a discussion 
in ref. |33|). At points where the right hand side of eq.(40) is regular it must 
coincide with the integral kernel of the analytic continuation of the semigroup 
exp(-Hj-). 



8 The semiclassical expansion 

A non-perturbative construction of QED is necessary for a demonstration of a 
classical limit of QED. This is so, because the conventional perturbative QED 
involves an expansion in the fine structure constant, where a — > oo when h — > 
( see a discussion of some problems concerning this limit in ref.Q). We could 
obtain the conventional perturbation expansion expanding eq.(27) in powers of 
the electric charge (see JjJ for a relativistic generalization of eq.(27)). In this 
section we obtain a quasiclassical form of the Feynman kernel (40) through an 
expansion in h (the semiclassical expansion of the Feynman kernel in quan- 
tum mechanics with a rigorous version of the Feynman integral is discussed in 
[[id] [^] (l9] ; for a version without the Feynman integral see ^3|). In such an 
approach the role of the two-point boundary value problem and its quantum 
correction (38) become visible. After a derivation of the quasiclassical form of 
the Feynman kernel other versions of the classical limit can be obtained through 
an expression of the time evolution by the Feynman kernel. In particular, one 
can see more explicitly the difference between QED and Boyer's stochastic elec- 
trodynamics 1 31]. In the derivation of the semiclassical expansion we follow our 
earlier paper [19]. The first step is a shift of variables 7 — * 7+C/3 and v — > v+zw, 
where v and w are real, whereas z and £ are complex numbers 
Lemma 8 

Let P(t) and w(r) be differentiable and their derivatives square integrable 
on the intervals [0, 1] and [0, t] respectively.Moreover, assume that /3 and w are 
equal to zero on the ends of the respective intervals. Then, with the assumptions 
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and notation of i) of Theorem 7 we have the following identity 



K = E 



ex p{^ ~ 3 So T,k,u C( fc > ^) 2 /3(t; k, v)C0(t; k, v)di 
- So £*,„ C(fc, «W(t; fc, i/) 7 (r; fc, i/)dr - £ £ r *r^Mv r (r) - ± £ £ r z^^r+ 
+*Er Jo ^(T,p r (r) + A£7- r « r w r (T);i,B + F,C)d(p r (r) + AfT r z r w r (r)) 

ex P^ _ £ Er, s e r e s Jo 0«(Pr( T ) + Acr r z. r w r (r),p s (r) + A(T r z r w r (T))dr 
-jE r Jo u r (p r (r) + Ao- r z r w r )dr^ 

(41) 

where C is defined in eq.(36) and 

A 5 (t, x; t, B + F, C) = ^(r, x; t, B, C) + T s {t, x) 

with 

•F,s(t,x) = ^2 y/4TTfi(k)cf(k)£(k, u)XakC(k, v)f5(h, u)g s (k, x) 

Proofuf z and C are real then eq.(41) (cp. the definition of 1Z in eq.(40)) is 
a version of the Cameron-Martin formula expressing the transformation of the 
Gaussian integral under a shift of variables; the validity of the formula (41) for 
complex z and C follows by analycity, this has been discussed first in ref. [Q 
and then in . 

Now, we choose z r = (Aoy) -1 and £(fc, v) = (A<7fc) -1 . Let us introduce a real 
field (see eq.(38) for a definition of A cl ) 

A*(t, x) = A c1 {t, x; t, B, C) + T{t, x) (42) 

and a real trajectory 

T £ — T 

x(r) = x — H — x + w(r) 

Let us choose wy in such a way that x r (r) is the solution of the Lorentz equation 
m r "37^ = ~e r V r E s e s (/)(x r , x s ) — V r u r (x r )— /• j o\ 

with the two-point boundary condition x(0) = x and x(i) = x'. For a sufficiently 
small time t (what we assume here) the two-point boundary value problem (43) 
has the unique solution (|4f|). 

The field A 1 is constructed from the amplitudes /3(r; k, v). Let us define 

£(r; k, v) = a(r; k, v) + /3(r; fc, v) (44) 
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where a is defined below eq. (34). Then, A 1 of eq.(42) is defined by £ 



A\r, x) = ^Hk)cf{k)£{k, v)Z( T ; k, v)g{k, x) 

Now, we choose /3(r; fc, f) in such a way that £(t; k, v) is the solution of the 
equation 



(dl + W (fc) 2 )£(r; fc, „) = ^y/(*) E fe^ 3 (fc, x(r)) (45) 

with the two-point boundary condition £(0; k, v) = X{k; v) and £(t; k, v) = 
Y(k;u). Eq.(45) implies the usual equation for a vector potential of a given 
particle current 

(l£ - A)^(x) = ^^S f (MrU) (46) 

r 

where 

«5 / (x,y) = l^/(fc) 2 3 (fc,x) 3 (fc,y) 

Eq.(46) is solved with the two-point boundary condition A 1 (t = 0, x) = B(x) 
and A 1 (t = t, x) = C(x). When /(fc) — > 1 then we obtain the classical elec- 
trodynamics of point charges. We wrote eqs.(43) and (46) without the 5- rcg- 
ularization. As a first step we insert again the <5-regularization and prove the 
scmiclassical expansion. Subsequently, we show that the 5-regularization can be 
removed. 

Theorem 9 

Let us assume that V is a cube and the number of modes is finite. Define 
the classical action (where x and £ solve the classical equations of motion (43) 
and (45)) 

Ss = Erlo (^(^) 2 - «r(xr(T)))dr- £ r , a /o e r e s Mx r ,x s )dT+ 

+ So (2m (T) 2 - §M*M*) 2 £(*, v?)dr + J2 r e T So Afc, x r (r))dx r (r) 

(47) 

then 

i) for a sufficiently small time t 

K s (t, C, x; 0, A, x!)exp{-ij-) = Q s (t, K) (48) 

is bounded in h 

ii) the limit Ss — > S as 8 — ► exists if t is small enough 
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iii) for a sufficiently small time £1$ can be expanded in a power series in the 
electric charge and the limit 5 — > of each term exists uniformly in fi. Hence, 
the limits 8 — > and h — ► in 0,5 (i, H) exist and are exchangeable. 

Proof: We expand K$ in powers of Vh according to eqs.(40)-(41).5' ( 5 of 
eq.(47) constitutes the singular part of the expansion. Equations of motion (43) 
and (45) result from the condition of the cancellation of the terms of order 
(y/K) -1 . There remains to bound the remainder. In the cube <7(fc,x) are the 
sinus and cosinus functions. In such a case the problem of the estimate of the 
remainder is very similar to the estimate of the remainder for trigonometric 
interactions done in ref. |]l9f (Theorem 6. 2). So, after a subtraction of the term Ss 
from Kg we expand the exponentials in cq.(41) in a power scries in the electric 
charge. The series is convergent for a sufficiently small time. We show that the 
singular terms of order (v^) _1 cancel identically in the series if the classical 
equations of motion (43) and (45) are fulfilled. The remaining terms in the 
series are already bounded in fi and 8. As a consequence, the series is convergent 
uniformly in Ti and 8. Hence, the limit 8 — > can be performed term by term in 
the convergent power series. The part ii) of the theorem is an easy consequence 
of the continuity of solutions of differential equations with respect to variation 
of parameters if the fields are regular and the time is small enough. The part 
iii) follows from our discussion of the expansion of Q$ in a power series. 

We can solve the linear equation (45) and insert the solution into eq.(43).In 
such a case we obtain an integro-differential equation describing particles inter- 
acting by Coulomb forces and additionally under the interaction of the retarded 
magnetic field produced by their movements. We would obtain this equation if 
we performed the classical limit in eq.(27) where the functional integration over 
the quantum electromagnetic field has been calculated. The double stochastic 
integral in the exponential of eq. (27) describes the quantum radiation reaction. 

Eqs. (43) and (45) lead to the classical expression for the radiation reaction 

[[46) j47fl • We have shown in Theorem 9 that the quantum radiation reaction has 
a classical limit for a sufficiently small time. It is well known that we encounter 
serious acausal effects with both the classical and quantum radiation reaction 
(see a clarification of some difficulties in ref. [[ID). However, these problems 
involving a large time behaviour are beyond the scope of this paper. 



9 Some approximations in QED 

It is quite difficult to obtain non-perturbative results in QED without further 
approximations. It is often argued (see e.g. p4|,[ 
in an atom is localized around the nucleus and the wave length involved in the 
light-atom interaction is large in comparison to the atom size then in quantum 
radiation theory we may neglect the spatial dependence of quantum electromag- 
netic field (such an assumption is also applied in mathematical physics papers 
[[u] jl5|]). Although this argument is convincing classically it would be very diffi- 



12 ) that because the electron 
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cult to justify it in a complete relativistic quantum theory of the electromagnetic 
field. In a relativistic theory the propagator D in eq.(27) is singular on the light 
cone. Hence, an expansion in coordinates of this propagator in eq.(27) could not 
be justified. We could consider the ultraviolet cutoff in D as an expression of 
the believe that the high frequency atom-photon interactions are not relevant 
to atomic physics. 

We intend to formulate this short range approximation in a particular case. Let 
us consider the formfactor 



We wish to compare the transition amplitude (27) with the one where the spa- 
tial dependence of D is neglected. Let us note that (see the propagator D of 
Proposition 6) 



D(x,x';|t-t'|)-£>(0,0;|t-t'|) = 

= Jo°° - x ') 2 (« + *c\t - r'l)- 2 ((es + iC \r r'|) 2 + (x - x') 2 ) 1 



Hence, this difference can be a bounded function of x — x'. 
Proposition 10 

Consider the transition amplitude (27) for a single particle in i? 3 . Assume 
that J s _2 <r(s)ds < oo.Then, for a sufficiently small time t 



where (, )o denotes the transition amplitude computed with the propagator 
D(0, 0; |r — r'|) where the spatial dependence is neglected. 

Proof: We note that under the assumptions of the Proposition the difference 
of the propagators D(Xah(r), Aerb(T'); |r — r'|) — D(0,0; |r — r'|) is a bounded 
function of b and r— r'. Then, in eq.(27) the exponential of the double stochastic 
integral J dbdb(D — Dq) (where the propagator without the spatial dependence 
is denoted Dq) can be expanded in a power series, which is convergent for a 
sufficiently small time. The first term J dbdb(D — D ) in the series (51) is 
estimated by the r.h.s. of eq.(51). Then, the estimate of the proposition is a 
consequence of the convergence of the series. 

The proposition says that under our assumptions an error of the neglect 
of the spatial dependence is of the order of the fine structure constant. Such 
a result can be considered as a first step of a more systematic investigation. 
The interesting result should hold true for arbitrarily large time. However, this 
would be quite difficult to prove. The difficulty results from the problem of in- 
tegrability for arbitrarily large time of the exponential of the double stochastic 
integral. The approximations to the double stochastic integral are interesting 
from the point of view of the quantum theory of radiation reaction. The well- 
known formula for classical radiation reaction of a point charge [flo] confirmed 




(49) 



(50) 



|($iX, U t <S>2X) - ($iX, U t $2X)o\ < const— t 



(51) 
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by approximate quantum calculations |^|[[30) can be obtained by a neglect of 
spatial dependence. In order to investigate the problem from the point of view 
of the Feynman integral (for a Hamiltonian formulation see [Q ) let us intro- 
duce the decomposition D = D° + D 1 for the zero temperature propagator and 
G = G° + G 1 for the finite temperature propagator (22), where the index zero 
means that the spatial coordinates are set equal to zero. In the infinite volume 
we have ( the factor | comes from the replacement of 8 tr by 5) 

D° (t) = -Sjt— / dkkf(k) 2 exp{-ickt) 



3 7T 



and 



2 Tic r°° f \ \ 
G°(t)ji = -Sji— / dkkf(k) 2 ( cth(-j3hck)cos{ckt) - isin(ckt) ) 

3 7T Jo V 2 J 

When f(k) — > 1 then ImD°(t)ji — > j^hSjiS'(t) (where Im denotes the imaginary 
part). Inserting this result into eq.(27) we obtain the standard expression for 
the radiation reaction (strictly speaking the integral (27) becomes divergent and 
requires a renormalization). We can take the radiation reaction into account 
explicitly by means of a transformation of paths. Let us introduce an operator 
ini 2 ([0,t]) 

p 2 

K = 1 + — — ImD° 

where by ImD° we denoted the integral operator defined by the kernel ImD . 
A transformation 

dh -► db = &dh 

expresses the functional integral (27) in terms of a new Brownian motion which 
contains the effect of the radiation reaction. Such a transformation is math- 
ematically correct if ImD is a Hilbert-Schmidt operator and moreover if /C 
is positive definite. For positive temperature it is useful to note that G — D is 
real. So, the heat bath does not contribute to the classical radiation reaction. In 
order to see the effect of the temperature we perform another transformation 
upon the functional integral. For any functional $ of the Brownian motion 
defined on the interval [0,t] we have the Cameron-Martin identity 



E[$(b)] = E 



for any f with a square integrable derivative. We choose now 

- = -A^(0,r) 

in eq.(24)( for simplicity restricted to a single particle). In such a case in eq.(24) 
$2(x + Aerb(t)) — * $2(<lt(x)) where q*(x) is the solution of the equation 

dq = -— A z (0, r)dr + Aadb(r) (52) 
m 
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Using eqs.(52) and (22) we can compute the correlation functions of q. In this 
way for a large temperature and f(k) = 1 we obtain 



^ ^1 
dr dr m 2 



plus less singular terms. So, the charged particle in equilibrium with a black 
body radiation behaves like a Brownian particle in a fluid of temperature T. This 
result is approximate, because we neglected the change of the covariance of the 
electromagnetic field caused by the Cameron-Martin transformation. We could 
approach this way in the framework of QED the model discussed first by Einstein 
and Hopf |l) and later on in refs.JHJ js^J.Lct us note that the behaviour of the 
correlation functions of the quantum particle in equilibrium with the radiation 
is in qualitative agreement with the requirements of the "quantum Brownian 
motion " due to Benguria and Kac Q( see also ref.Q). 



10 Discussion and outlook 

In this paper we concentrated mainly on the precise formulation of regularized 
QED dynamics. In some cases we were able to show that the expressions were 
finite only for a small time. Such a requirement resulted from a demand of a 
convergence of a perturbation series. In principle, we can prolong the dynamics 
to a larger time using the group property of the time evolution. In practice, 
the multiple integration involved in such a procedure may be inconvenient. The 
results of sec. 7 may indicate a way out of this difficulty. The formula for the 
propagation kernel of a harmonic oscillator holds true for arbitrarily large time 
(except of the moments equal to the classical periods) in spite of the fact that 
the functional integral looses its sense already when t tends to the first pe- 
riod. In another approach to the harmonic oscillator we have obtained in sec. 3 
the dynamics for an arbitrarily large time through an application of stochastic 
differential equations (see eq.(5)). We have discussed the method of stochas- 
tic differential equations for non- linear quantum systems in refs. fig] |5^] . This 
method could be applied to QED. The first step has been done in eq.(52). A 
non- linear transformation (discussed for quantum mechanics in |pH ) can be ap- 
plied to formulate exactly the dynamics in QED as a random dynamics. We 
intend to discuss this method in a subsequent publication. 

The dependence on the formfactor (ultraviolet cutoff) constitutes another 
unsatisfactory aspect of the theory. A non-perturbative removal of the cutoff 
is a difficult problem. However, some results may be insensitive to the cutoff. 
One may believe that phenomena on a large time and space scale involving low 
energies should not depend on the ultraviolet cutoff. We consider the formula- 
tion of this paper as a substantial step towards a control over approximations 
made in QED, quantum optics and quantum radiation theory. In principle, in 
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sec. 9 we were able to estimate an error involved in a dipole approximation. In 
sec. 8 we are able to estimate subsequent terms in the semiclassical expansion. 
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